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Abstract: The Pierre Auger Observatory is exploring the potential of the radio detection tech- 
nique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger En- 
gineering Radio Array (AERA) addresses both technological and scientific aspects of the radio 
technique. A first phase of AERA has been operating since September 2010 with detector stations 
observing radio signals at frequencies between 30 and 80 MHz. In this paper we present com- 
parative studies to identify and optimize the antenna design for the final configuration of AERA 
consisting of 160 individual radio detector stations. The transient nature of the air shower signal re- 
quires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not 
widely discussed in antenna literature, we review the relevant antenna characteristics and enhance 
theoretical considerations towards the impulse response of antennas including polarization effects 
and multiple signal reflections. On the basis of the vector effective length we study the transient 
response characteristics of three candidate antennas in the time domain. Observing the variation of 
the continuous galactic background intensity we rank the antennas with respect to the noise level 
added to the galactic signal. 

KEYWORDS: Cosmic Rays; Air Showers; Radio Detection; Antennas; AERA; Pierre Auger 
Observatory. 
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1. Introduction 

The Pierre Auger Observatory [1] is a hybrid detector for the observation of cosmic rays above 
~ 10 18 eV. The atmosphere of the Earth is used as a calorimeter measuring the particle shower that 
evolves after the penetration of a primary cosmic ray. These air showers are observed with an array 
of 1660 ground based particle detectors covering an area of 3000 km 2 . The array is overlooked by 
27 optical telescopes [2] which are sensitive to the fluorescence light emitted by nitrogen molecules 
which were excited by the charged particles of the passing air shower. The combination of both 
detection techniques allows for a precise determination of the energy and arrival direction of cosmic 
rays and gives information on the chemical composition of the cosmic ray flux [3,4]. 

Besides the established detection techniques, the Pierre Auger Collaboration is exploring the 
possibility of detecting extensive air showers via radio pulses that are generated as the showers 
develop in the atmosphere [5-7]. The radio signal strength promises a quadratic scaling with the 
energy of the cosmic ray, high angular resolution in the reconstruction of the air shower axis, and 
sensitivity to the nature of the primary particle. In combination with an almost 100% duty cycle 
these attributes make a radio system a candidate for the next generation of ground-based air shower 
detectors. 

The emission of electromagnetic radiation from air showers in the MHz frequency regime 
was first observed by Jelley and co-workers in 1965 [8]. It was found that air showers emit an 
electromagnetic pulse in the direction of the shower propagation. The observation of the wavefront 
with an array of individual antennas at different positions with respect to the shower axis should 
allow a reconstruction of the properties of the air shower and the corresponding primary cosmic ray. 
In the following years, progress was made with experiments reporting air-shower observations in 
a frequency range from 2 to 550 MHz [9, 10]. The realization of a comprehensive radio detector, 
however, was not feasible until the appearance of fast digital oscilloscopes in the past two decades. 
Since 2005 the CODALEMA [11] and LOPES [12] experiments have succeed in detecting air 
showers up to energies of 10 18 eV. They confirmed the predicted importance of the Earth's magnetic 
field for the generation process of the radio pulse [13, 14]. 

The Auger Engineering Radio Array (AERA) is a radio detector situated at the Pierre Auger 
Observatory. AERA will instrument a sensitive area of 20 km 2 with 160 detector stations and is thus 
the first detector with sufficient collecting area to make possible the measurement of radio signals of 
air showers beyond 10 18 eV. The layout of AERA shown in Fig. |l| It features a varying spacing 
between the detector stations which is intended to maximize the number of recorded events over 
a wide energy range from roughly 10 17 to 10 19 eV at a rate of several thousand air showers per 
year [15]. AERA is co-located with Auger fluorescence telescopes and its surface detector. Hence, 
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Figure 1 : Layout of the Auger Engineering Radio Array with its radio detector stations embedded 
in the grid of particle detectors in the north-western part of the Pierre Auger Observatory. AERA 
has an extent of ~ 20 km 2 and consists of 160 radio detector stations that will be deployed in three 
stages. The first stage of AERA is operating since September 2010. The Central Radio Station hosts 
the central data acquisition of AERA. 

the Pierre Auger Observatory offers the unique possibility to study radio emission even at large 
distances from the shower axis. The first stage of AERA consists of 21 autonomous detector stations 
forming the dense core seen in Fig. [j] and has been operating since September 2010. 

AERA is a self-triggered radio array: at each detector station the radio signal is observed 
continuously and a trigger decision is formed to select air showers which produce a transient signal 
of a few times 10 ns in length. Individual trigger decisions are collected in a data acquisition system 
where signal patterns are investigated that match air showers recorded in multiple detector stations. 

The radio detector stations provide ultra-wideband reception of radio signals from 30 to 80 
MHz, a bandwidth locally free of AM and FM band transmitters. The continuous radio signal is 
sampled at a rate of 1 80 MHz and processed with a field-programmable gate array (FPGA) to form 
trigger decisions. A review of the readout electronics is given in Ref. [16]. Fig. shows a picture of 
a radio detector station. 

The data recorded by the detector stations is a convolution of the radio signal and the response 
of the readout electronics. To create a measurement that is independent of the detector setup, the 
impact of the individual hardware components has to be removed from the AERA data. Here, the 
antenna deserves special attention as its frequency response is highly non-linear and depends on 
the incoming direction and polarization of the recorded signal. Hence, a precise knowledge and 
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Figure 2: A radio detector station of the first stage of AERA. A logarithmic -periodic dipole antenna, 
called Small Black Spider, is used as sensor for the radio emission of extensive air showers. The 
antenna is read out with digital oscilloscopes and a trigger decision on air shower signals is performed 
by an FPGA. A photovoltaic power system enables autonomous operation of the detector station and 
supplies the readout electronics, the GPS system used for timing, and the communication system. 



proper use of the relevant antenna characteristics is required to obtain a calibrated measurement of 
the radio emission from air showers. 

The antenna and its amplifier determine the signal to noise level that is obtained in air shower 
observations. An optimal noise performance of the antenna will maximize the sensitivity of the 
radio detector stations to air shower signals and the efficiency of the detector to cosmic rays. 

In this article we present our studies to characterize and evaluate candidate antennas for the 
next stage of AERA. The article is structured as follows: in section 2 the candidate antennas are 
presented. In the third section we focus on the antenna theory needed for understanding the reception 
of transient signals. Having identified the relevant antenna characteristics we show in section 4 
how to access them in simulation and measurement. In section 5 we compare the response of the 
tested antennas to transient signals on the basis of antenna simulations. Preferably, an antenna for 
the detection of air showers should introduce only minor distortions to the recorded signal shape. 
The successful observation of air showers will be determined by the signal to noise ratio obtained 
in measurements. Therefore we present in section 6 comparative measurements of the variation 
of the galactic noise level performed at the Nancay Radio Observatory, France, which allow us to 
discriminate the candidate antennas with respect to their noise performance. 
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2. Antennas for the Detection of Radio Signals from Cosmic-Ray Induced Air Show- 
ers 

In this section we present three antenna models that have been evaluated for the setup at AERA. 
The three investigated models address the task of radio detection of air showers with different 
concepts. All antennas presented are the result of several stages of development, taking into account 
the experience gained in smaller radio detection setups both in Europe and at the site of the Pierre 
Auger Observatory. This is necessary as the environmental conditions of the Argentinian Pampa 
impose special demands on the antenna structure for instance through wind loads of up to 1 60 km/h. 
Hence, the durability and consequent costs for maintenance directly impact on the success of an 
antenna model, especially for a radio detector design with a large number of detector stations. 

With respect to the electrical properties of the radio sensors it is useful to consider the antenna 
as an integral combination of the metallic structure capturing the signal and the first low noise 
amplifier (LNA). Whereas the structure determines the directional properties of the antenna, the 
ultra-wideband reception of the antenna is ruled by the combination of both elements. 

Requirements on the directional properties of the antenna are imposed by the widespread 
layout of the radio detector array. At each detector station the full sky needs to be observed so 
omni-directional antennas are used. To measure the polarization of the radio signal the placement of 
at least two perpendicular antennas is required at each detector station. 

Currently, the properties of the radio pulse and its generation mechanism are subject to research 
beyond its capabilities as a tool to detect the cosmic ray. Hence the antennas have to be sensitive in 
a broad frequency range to allow for a maximum detail of the observation. At AERA the bandwidth 
is limited by the presence of AM band transmitters below 30 MHz and FM band transmitters mainly 
above 80 MHz. 

2.1 The Small Black Spider Antenna 

Logarithmic periodic dipole antennas (LPDAs) are being used for the first stage of AERA. The 
logarithmic periodic principle assembles a series of half wave dipoles of increasing length to keep the 
radiation resistance of the antenna constant over a wide frequency range. LPDAs were first adapted 
to the needs of radio detection for the LOPES-STAR experiment [17]. The 'Small Black Spider' 
realizes the LPDA principle as a wire antenna and is shown in Fig. 0. The antenna integrates two 
independent antenna planes in the same mechanical structure which has a dimension of 4 x 4 x 3.5 
m 3 and a weight of 18 kg. For transportation purposes the design of the Small Black Spider includes 
a folding mechanism in the antenna structure. This allows one to assemble the antenna completely 
in the laboratory and make it operational within 15 minutes at the detector site. 

The lengths of the shortest and longest dipoles determine the available frequency range and 
have been adapted to the AERA frequency and for the Small Black Spider. A slightly enhanced 
sensitivity of the antenna to the top direction is obtained as the amplitudes fed from the individual 
dipoles into a common wave guide add up constructively at the footpoint where the antenna is read 
out. The footpoint of the Small Black Spider is at the top of the structure which is in a height 
of ~ 4.5 m when installed in the field. Although the footpoint is the optimum position for the 
amplifier such placement is not feasible due to maintenance constraints. Instead a matched feed 
of the footpoint into a 50 £2 coaxial cable is obtained using a transmission line transformer with 
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Figure 3: The Salla antenna during test measurements at the Nancay Radio Observatory. The Salla 
realizes a Beverage antenna as a dipole loop of 1.2 m diameter. The antenna is read out at the top 
with an LNA. The amplified signal is guided through the antenna structure to the bottom. 

a 4:1 impedance ratio. This is possible because of a constant antenna impedance of 20011 within 
the frequency range. The coaxial cable guides the signal to the LNA at the bottom of the antenna. 
Here, the LNA includes filter elements at its input to further ensure the frequency selectivity of the 
antenna to the AERA band. Details of the development of the LNA are presented in Ref. [18]. 

For the first stage of AERA, 30 Small Black Spiders have been produced including spares 
and antennas for test setups. Within one year of placement in the field the antennas have proven 
to be both robust against the environmental conditions of the Argentinian Pampa and sensitive to 
radio signals from cosmic ray induced air showers. A detailed description of the antenna is given in 
Ref. [19]. 

2.2 The Salla Antenna 

The short aperiodic loaded loop antenna 'Salla' realizes a Beverage antenna [20] as a dipole loop of 
1.2 m diameter. The Salla has been developed to provide a minimal design that fulfills the need for 
both, ultra-wideband sensitivity, and low costs for production and maintenance of the antenna in a 
large scale radio detector. The compact structure of the Salla makes the antenna robust and easy to 
manufacture. A picture of a Salla showing its two polarization planes is displayed in Fig. 0. 

Beverage antennas include a resistor load within the antenna structure to give a specific shape 
to the directivity. In the case of the Salla a resistance of 500 Q. connects the ends of the dipole arms 
at the bottom of the antenna. The antenna is read out at the top which is also the position of the 
LNA. While signals coming from above will induce a current directly at the input of the amplifier, 
the reception from directions below the antenna is strongly suppressed as the captured power is 
primarily consumed within the ohmic resistor rather than amplified by the LNA. The resulting 
strong suppression of sensitivity towards the ground reduces the dependence of the antenna on 
environmental conditions which might vary as a function of time and are thus a source of systematic 
uncertainty. With the inclusion of an ohmic resistor the Salla especially challenges its amplifier as 
only ~ 10% of the captured signal intensity is available at the input of the LNA. Proper matching 
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Figure 4: A Butterfly antenna installed at the site of the Nancay Radio Observatory. The antenna 
construction aims at a close integration with a box containing the readout electronics and a center 
pole which is also used to carry an antenna for wireless communication. 

between the antenna structure and the LNA is realized with a 3:1 transmission line transformer. The 
structure of the Salla creates a sensitivity which is flat as a function of frequency. 

Salla antennas are used at the Tunka radio detector [21]. They have been tested in a radio 
detection setup at the Pierre Auger Observatory and within the LOPES experiment. A detailed 
description of the Salla antenna is given in Ref. [22]. 

2.3 The Butterfly Antenna 

The 'Butterfly' is an active bowtie antenna and the successor of the active short dipole antenna [23] 
used for the CODALEMA radio detector. The Butterfly has dimensions of 2 x 2 x 1 m 3 and is 
constructed with fat dipoles. The dipole signals are fed directly into the two input channels of 
an LNA at the center of the antenna which constitutes the active antenna concept. The dipoles 
are self-supporting and their hollow construction reduces the sensitivity to heavy wind loads. The 
Butterfly antenna is shown in Fig. |]. 

The Butterfly explicitly uses the presence of the ground to enhance its signal. The center 
of the antenna is installed at a height of 1.5 m. Here, the direct wave and the wave reflected 
on the ground add constructively in the antenna throughout most combinations of frequency and 
incoming directions. Ultra-wideband sensitivity is obtained by designing the input impedance of 
the LNA to depend on the impedance of the fat dipole structure as a function of frequency. In this 
way, the sensitivity of the dipole has been optimized allowing for an efficient detection also of 
wavelengths that are much longer than the dimension of the antenna structure. The LNA itself is a 
application-specific integrated circuit (ASIC) and does not require a transformer since its input is 
differential [23]. 

In the current extension of the CODALEMA experiment [24] 33 Butterfly antennas are deployed. 
The antenna was used successfully to observe cosmic rays in one of the pioneering setups at the 
Pierre Auger Observatory [25]. Details of the Butterfly antenna are given in Ref. [26]. 



Figure 5: The spherical coordinate system with the antenna structure in the center. The origin of 
coordinate system is located in the XY-plane below the antenna which takes the placement of the 
antenna above a ground plane into account. Depicted is a logarithmic periodic antenna structure. 
The zenith angle 6 is counted from the top, the azimuth angle counterclockwise from the x-axis 
of the coordinate system. A specific direction (6, 0) is considered as the incoming direction of a 
signal. The vector of the effective antenna length H for the specified direction is given. H as well as 
the vector of the electric field (not depicted) are contained in the plane spanned by the unity vectors 
eg and e$. 

3. Antenna Theory 

To perform a calibrated measurement of the radio emission from cosmic ray induced air showers the 
impact of the detector and especially of the antenna needs to be unfolded from the recorded signals. 

The goal of this section is to introduce theoretical aspects of ultra-wideband antennas needed to 
describe the interrelation between the measured voltage V(t) responding to an incident electric field 
E(t). Here, we aim to unify the calculations for the diverse antennas described in Sec. |[ Having 
identified the relevant quantities, a closer evaluation of the antenna models will take place in the 
following sections of this article. 

3.1 The Vector Effective Length 

For antenna calculations it is convenient to choose a spherical coordinate system with the antenna 
in its center as depicted in Fig. |5|. In this coordinate system the electric field of a plane wave that 
arrives from a given direction (0,0) at the antenna will be contained in the plane spanned by the 
unity vectors eg and only. The electric field can be written as a two-component vector and is 
called the instantaneous electric field. Its two components denote two independent polarization 
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directions which vary as a function of time: 



E(t) = e e E g (t) + e^E^(t) . (3.1) 

The mapping between the voltage response V(t) induced over the antenna terminals and the 
electric field E{t) is represented by the vector effective length (VEL) H{t) of the antenna [27]. As 
displayed in Fig. |5] the VEL is a two-component vector in the antenna-based coordinate system as 
well: 

H{t)=e e H (t)+e$H$(t) . (3.2) 

Here, Hq (?) encodes the response characteristics of the antenna to the component of the incident 
field in eg -direction and (?) accordingly. The VEL contains the full information on the response 
of an arbitrary antenna structure to an arbitrary plane wave signal. 

The response voltage of the antenna to a single polarization direction of the electric field is 
obtained by the convolution of the field and the VEL. For instance in the -direction the response 
is calculated as: 

V$ (t)=Ht(t)*Et(t) , (3.3) 

where the symbol V marks the convolution transform. The total antenna response is obtained as 
superposition of the partial response voltages V§ and Vg to the two independent polarizations of the 
electric field [28]. Using a vectorial form we can write conveniently: 

V(t)=H(t)*E(t) . (3.4) 

Up to now we have treated the antenna response calculation in the time domain. However, the 
antenna characteristics contained in H(t) are usually accessed in the frequency domain rather than 
the time domain. We define the vectorial Fourier transforms of the quantities, e.g.: 

U{a) = e e T(H e {t))+e$F(H$(t)) . (3.5) 

and V, £ accordingly. The script letters indicate complex functions of the angular frequency 
CO = 2nv belonging to the frequency v. 

The convolution theorem allows the convolution of functions to be performed as a point wise 
multiplication of their Fourier transforms. Hence, the voltage response of Eq. ^3] can be treated in 
the frequency domain as follows: 

V((0) =H((0) -£{a) . (3.6) 
The voltage response in the time domain follows from the inverse Fourier transform V{t) = 



F 1 (V(ft))). It should be noted that Eq. represents a condensed way to calculate the antenna 
response to arbitrary waveforms. 

3.2 Polarization 

In Fig. |5] the antenna-based spherical coordinate system along with an exemplary antenna structure 
is shown. The pictured vector H can either be understood as the VEL at a certain point in time 
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Figure 6: The VEL in the case of the depicted logarithmic periodic dipole structure derived from 
simulations with the numerical antenna simulation program NEC-2 [29] at 75 MHz. For a single 
frequency the VEL can be expressed as a vector field in the spherical coordinate system of the 
antenna. 

H(t = ?,-), or at a given frequency % (oo = COj). In the latter case the components of % are complex 
functions: 

% = ee%e+?0%0 , (3.7) 
= 3e\He\j* H >+2*\Ut\e Uhl * . (3.8) 

At a specific frequency the VEL in Eq. [£7] resembles a Jones vector [30] which is commonly used 
to describe the polarization of light. Of special interest is the phase difference between the two 
components: 

A<?>H = <PH^ - <PH e ■ (3-9) 
Using this phase difference we separate Eq. into a global and a relative phase: 

U = e'vxe (e e \H 9 \+e$ |7^|e ! ' A<p «) . (3.10) 

If the phase difference A(f>u at a given frequency CO is a multiple of % 

Aq>H = nn, n = ...,-1,0,1,... (3.11) 

the maximum sensitivity of the antenna is reached for the reception of a linear polarized signal. In 
this case an intuitive picture of the antenna can be drawn. Following the Jones calculus, we rewrite 
Eq. |3.10| omitting the global phase: 



H' = \H\ (eg cos a + e<j, sin a) 



(3.12) 



where 

a = (-l) n arctan(^4) (3.13) 
\rLe\ 

denotes the angle of the VEL axis in the eg-e^ -plane counted counterclockwise from eg. The VEL 
%' is pictured in Fig. || for a single direction on the unit sphere. In Fig. ^ %' is displayed for a set of 
incoming directions. 

The displayed characteristics were accessed via simulations for the displayed logarithmic- 
periodic antenna structure. The length of the vector changes as a function of zenith angle which 
denotes the directionality of the antenna. The LPDA is most sensitive to the vertical (zenith) 
direction. 

In the setup shown in Fig. the VEL vanishes when approaching the x-axis. For the given 
antenna structure, the electric field of a wave incoming along the x-axis will have no components 
parallel to the dipoles of the antenna and cannot be detected. 

For a simple antenna structure as in Fig. the VEL is aligned with the projection of the 
antenna dipole on the unit sphere for a given direction. For any incoming direction of the wave 
a configuration of the electric field and the effective antenna length exists such that no signal is 
detected. This is referred to as polarization mismatch. 

For single frequencies, the VEL is thus a vector field of 2-dimensional complex vectors in the 
antenna-based spherical coordinate system. The VEL thus has three major dependencies: 

n((o) = n((o,d,(t>) . (3.i4) 

In the general case of an elliptical polarization where A<p^ ^ n% the VEL can be presented as 
an ellipse on the unit sphere rather than a vector. However, in the case of the simulated antenna, the 
VEL is a vector. The omission of a global phase in Eq. [3.12| enables the inspection of the polarization 
of an antenna at single frequencies. The wideband characteristics of antennas are contained in the 
development of % as a function of frequency. 

The VEL is related to the antenna gain and the antenna directivity. In Appendix [A| we briefly 
discuss its relationship to these more commonly used quantities. 

In the case of the antenna models discussed in Sec. @, two rotated antennas are assembled with 
the same hardware structure. In Appendix [B] we show how the VEL can be applied to reconstruct 
the 3-dimensional electric field vector of the signal recorded in such antenna setups. In Appendix [B] 
we also relate the VEL to the Jones antenna matrix used in radio polarimetry [31]. 

3.3 The Realized Vector Effective Length in a Measurement Setup 

The VEL as defined in Sec. f57T] relates the incident electric field to an open circuit voltage at the 
antenna terminals which we will refer to as V c in the following. In an actual measurement setup 
the antenna will be read out at a load impedance. In this section we will focus on the impact of the 
readout system on the measured signal. 

3.3.1 The Thevenin Equivalent Antenna 

In Fig. [7| (left) the Thevenin equivalent circuit diagram for a simple measurement situation is 
displayed. The antenna is read out introducing a load impedance Z/, in addition to the antenna 
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Figure 7: Left: Thevenin equivalent of an active antenna used for reception. Right: Thevenin 
equivalent of an antenna used for reception with intermediate transmission line and transmission 
line transformer. Please refer to the text for a description of the symbols. 



impedance Z4. The voltage measured over the load impedance follows from the voltage divider 
relation: 

V L = -^—V oc =pVoc ■ (3.15) 

The impedances Za and are complex functions of the frequency ft). The situation of conjugate 
matching (cf. Appendix |A|) is obtained when the readout impedance is the complex conjugate of the 
antenna impedance Zl = Z* a . 

If the complex transfer function p in Eq. |3.15| is included in the formulation of the measurement 



equation (Eq. |3.6| ), the VEL is referred to as the realized or normalized VEL % r [32]: 

V L = pn-£ = H r -£ . (3.16) 

The transfer factor p as given in Eq. [3.15| is sufficient to calculate the response voltage in the case of 
the Butterfly antenna (cf. Sec. 0). For the other antennas, more elaborate diagrams are required. 

3.3.2 Transformers for Impedance Matching 

The transfer function p may include more complex setups than the one displayed in Fig. [7J (left). 
The small aperiodic loop antenna Salla includes a transformer to realize a better matching between 
the antenna impedance Za and the impedance of the readout amplifier Z^. 

We denote the impedance transformation ratio of the transformer by r. In the circuit diagram 
shown in Fig. [7] (left) the load impedance due to the combination of readout impedance and 
transformer is then r • Zi . Along with impedance transformation the transformer changes the voltage 
that is delivered to the readout impedance. Here an additional factor 1 / yfr needs to be added [33]. 
Hence, if an ideal transformer is used to optimize matching, the transfer factor p becomes: 

y/r Z A + rZ L 

3.3.3 Intermediate Transmission Lines 

The calculations performed in Sec. |3.3. 1| and |3.3.2| implicitly assumed that electric distances between 
the position of the impedances are short in comparison to the wavelength processed. However, the 
Small Black Spider LPDA uses a transformer to feed the antenna signals into a coaxial cable which 
guides them to the first amplifier. With increasing length of the coaxial cable, the direct current 
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Figure 8: The setup for calibration measurements of the Small Black Spider LPDA at the site of 
AERA. Left: Sketch of the measurement setup. A calibrated transmitter antenna is moved on a 
circle around the antenna under test to access different zenith angles. Only the read out plane of the 
tested antenna is displayed. Middle: A picture of the setup. In the lower left the AERA station used 
for the calibration measurement is visible. At a distance of ~30 m a balloon carries the calibrated 
transmitter antenna. The position of the transmitter is fixed by a three-legged rope system which 
provides the movement on the circle and a parallel orientation of the transmitter and receiver antenna. 
Right: Picture of the calibrated biconical antenna used as transmitter. The biconus has a length of 
1 .94 m. The antenna is carried by a balloon filled with 5 m 3 helium and is fed by a coaxial cable 
running to the ground. The overall weight lifted by the balloon is ~4.5 kg. 



approximation becomes invalid and propagation effects need to be taken into account. A circuit 
diagram for this setup is displayed in Fig. [7] (right). In Appendix we derive a transfer factor p that 
enables calculations of multiple signal reflections between readout impedance and antenna to be 
made as a single step in the frequency domain. It is given by: 

- (1+rJe "' (3.18) 



r Z A + rZ„ el*>'-T A T L 

Here, Z t i is the characteristic impedance of the intermediate transmission line and I' its electrical 
length. The complex propagation constant y per unit electrical length includes the attenuation loss 
along the transmission line. Tl and are the voltage reflection coefficients from the transmission 
line to the load and from the transmission line to the antenna respectively. 

The results for the less complex setups discussed in the previous subsections are included in Eq. 
|3.18| , e.g. when the transmission line is short: V — > or the transformer ratio is r = 1. Hence, Eq. 
13.181 unifies the calculation for the three antennas discussed in Sec. 0. 

In the case of the Small Black Spider LPDA, the electrical line length is /' 4.4m. With respect 
to the discussion in this section special care was taken during the design of the amplifier impedance 
to match the 50 Q. transmission line. Hence, Eq. [3 .18| introduces only slight changes to the signal 
shape and a time delay due to the length of the transmission line. 
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3.4 Signal Amplification 

With Eq. |3.16| the incident electrical field is related to the voltage Vl over the readout impedance of 
the antenna. For the antennas being considered this readout impedance is the input impedance of an 
LNA. Amplifiers are characterized by the complex scattering or S-parameters [34]. The S-parameter 
S21 is the ratio of the amplified to the incoming voltage amplitude. It should be noted that the 
amplification characteristics expressed by S21 implicitly assume that the amplifier is operated within 
a system with the same impedance as during the S-parameter measurement. Looking e.g. at Fig. [7] 
(left) this assumption is not fulfilled: the antenna impedance Za replaces the generator impedance 
used during the S-parameter measurement and may itself be frequency dependent. In Appendix [D| 
we show that, in this case, a corrected amplification factor S21' — also referred to as voltage gain — 
can be used. This is normalized to the voltage over the amplifier input impedance Vl during the 
S-parameter measurement: 

S21 

S21< = — . (3.19) 

Here, S 1 1 corresponds to the voltage reflection coefficient at the input of the amplifier during the 
S-parameter measurement. In the case of the antennas discussed, S21' relates the voltage V a at the 
amplifier to the voltage amplitude fed into a well defined 50 £2 system: 

V a = S2l'-V L . (3.20) 



This amplified voltage is used for further signal processing. V a is of special interest as signal to noise 
ratios are essentially fixed in a readout chain after the first amplification. The noise performance is 
not discussed here but will be determined experimentally in Sec. ^for the discussed antennas. 

Especially for antennas where the LNA is directly integrated into the antenna structure, e.g. for 
the Butterfly antenna, it is useful to define a VEL that includes the amplification of the signal: 

n a = s2\'-n r . (3.2i) 

1-L a is convenient to access in antenna calibration measurements as we will show in Sec. |[ Using 
Eqs. |3.20| and |3.18| allows us to take the impact of the amplifier and the readout system into account 
when studying antenna characteristics simulated with the NEC-2 program in the next sections. 



4. Measurement of the Characteristics of the Small Black Spider LPDA 

In Sec. we have recognized that the vector effective length is the central antenna characteristic 
needed to perform calibrated measurements of the electric field strength of incoming radio signals. 
The VEL can be calculated in simulations. However, it is desirable to perform the unfolding of the 
impact of the detector on the recorded signal on the basis of measured antenna properties. 

In this section we present calibration measurements we have performed to access the VEL of a 
Small Black Spider LPDA installed as part of the first stage of AERA. In our measurements we focus 
on the zenith angle dependence of the VEL as a function of frequency. In preceding measurements 
we have confirmed that the azimuthal dependency of the VEL follows a simple sinusoidal function 
as expected for dipole-like antennas such as the Small Black Spider LPDA [35]. 
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4.1 Antenna Calibration Setup 

So far, directional properties of antennas for the detection of cosmic rays have been measured with 
down-scaled versions of the antenna under test [36] or have been constrained to a measurement of 
the scalar amplitude transfer without the vectorial phase information [37]. 

The experimental challenge in a calibration measurement with a full scale radio antenna is to 
fix a distance between transmitter and receiver antenna such that the wave emitted by the transmitter 
sufficiently approximates the plane wave condition at the tested antenna. It should be noted that the 
usual approximation at which distance rff this far field condition is fulfilled is given by rff > 2 8 2 /X, 
where 8 the largest dimension of the transmitting antenna. This relation is valid only if 8 is larger 
than the wavelength X [38]. This is typically the case for dish antennas but not for the radio antennas 
discussed in this article. 

With respect to the longest wavelength of the AERA bandwidth (10 m), distances between 
transmitter and the tested antenna of 3 A are realized in our calibration. Following the discussion 
in Ref. [39] we can estimate the impact of near-field components still present in the calibration 
measurement. At a distance of 3A these will cause a variation of the power angular distribution of at 
most ±0.5 dB when compared to a measurement performed at much larger distances. 

In Fig. |H| an overview of the calibration setup is given. To access large distances at small zenith 
angles above the antenna under test, a balloon is used to lift a calibrated transmitter antenna [40]. 
Ropes constrain the movement of the transmitter to a circle around the tested antenna and ensure 
a parallel orientation of the two antennas. In this way the -component of the VEL can be 
accessed for various zenith angles. 

The transmission measurement from the biconical antenna to the Small Black Spider is per- 
formed using a vector network analyzer [41]. The network analyzer simultaneously feeds the 
transmitter antenna and reads out the tested antenna. The signal delivered by the vector analyzer 
is adjusted to appear > 30 dB above the ambient radio background recorded by the tested antenna 
throughout the measurement bandwidth from 30 to 80 MHz. The impact of the coaxial cables needed 
for the connections are removed from the data by including them in the null calibration of the vector 
network analyzer prior to the measurement. The amplifier of the Small Black Spider is included 
in the transmission measurement. Hence, the setup allows us to measure % a ,$ (fi), $ , = 270°) as 
discussed in Eq. [3.2 1| of Sec. ^4]. 

4.2 Simulated Calibration Setup 

To cross-check the calibration measurement procedure we performed simulations of an equivalent 
setup using the numerical antenna simulation tool NEC-2 [29]. The simulated calibration setup 
includes a model of the Small Black Spider LPDA as well as a model of the biconical antenna used 
as a transmitter. Both antenna models are placed in the simulation according to the geometries 
existing in the actual calibration setup. The sketch in Fig. || (left) is generated from a NEC-2 
simulation of a transmission measurement. 

The simulation model of the transmitter antenna is excited by placing a voltage source at its 
footpoint. The NEC-2 simulation then calculates the power emitted by the transmitter and the 
consequent open terminal voltage V c induced in the structure of the Small Black Spider. Following 
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the discussion in Sec. ^3] we process the terminal voltage to give the power delivered into a 
Z t i = 50 £2 system which corresponds to the coaxial cables connected to the network analyzer. 

In the field measurements, not all power that is delivered by the signal source of the network 
analyzer is accepted by the transmitter due to impedance mismatch. This has to be accounted for 
in the simulation by dividing the radiated power calculated by NEC-2 with the measured power 
acceptance of the transmitter antenna. 

In the measurements at the AERA site, signal reflections from the ground are included. We 
thus use the option of NEC-2 to model a ground plane in the simulation using the Fresnel reflection 
coefficients. The reflection coefficients depend on the relative permittivity £ r and the conductivity 
a of the soil. With a = 0.0014 n _1 m _1 we assume a low conductivity which has been confirmed 
in initial test measurements at the AERA site. For low conductivities we find £,. = 5.5 to be a 
reasonable choice for the relative permittivity in typical ground scenarios [42]. The resulting ground 
exhibits a relatively low reflectivity. 

As in real measurements, we use the simulation to yield the ratio of signal received by the 
tested antenna to the signal used to operate the transmitter. Simulated and measured data are then 
processed equally in further analysis. 

4.3 Transmission Equation and Data Processing 

For each zenith angle accessed in the setup we measure the S-Parameter S21 as a function of 
frequency. S21 is the complex ratio of the voltage amplitude V a delivered by the tested antenna and 
the amplitude V g delivered from the signal generator to the transmitter: 

S21 = -f . (4.1) 

The voltage V a is the response of the tested antenna to the electric field £^ caused by the transmitter 
antenna: 

V a = U a j£\ ■ (4.2) 

Due to the configuration of the setup, the electric field is contained in the e§ -direction of the antenna- 
based coordinate system. In antenna theory (e.g. Ref. [38]) the electric field of a transmitting 
antenna at a distance R is given. For our calibration setup it is: 

Z\ = -iZ*^fi%e- ie > R l c , (4.3) 

where is the VEL, 2q denotes the current in the transmitter antenna as depicted in Fig. [7], and c 
is the speed of light. When Eqs. and |3] are inserted in Eq. |TT] we obtain a complex form of the 
Friis transmission equation: 

S2\ = -iZ ^-^-U\n a ^e- imR ' c . (4.4) 

To access the desired VEL % a ,s of the antenna under test, the characteristics of the calibrated 
transmitting antenna %\ have to be applied. In our case these are given in terms of realized gain G r . 
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Figure 9: The amplified vector effective length of the Small Black Spider LPDA for the zenith 
direction as a function of frequency in measurement and simulation. The uncertainty of the 
measurement is indicated by the shaded area and is dominated by the systematic uncertainties 
of the calibration of the transmitting antenna of 0.7 dB and by the precision of the transmission 
measurement, here 0.6 dB. 

The realized gain refers to the transfer of signal power rather than signal amplitude and includes the 
reflections at the input of the transmitter antenna. In Appendix |E| we derive that: 



The transmitter calibration G' ^((o) is given by the manufacturer of the antenna as a function of 
real numbers. Such simplification is acceptable if the transmitter antenna introduces only minor 
distortions to the phasing of the signal within the measurement bandwidth. This has been verified in 
preceding test measurements for the biconical antenna [35]. 

The combination of Eqs. and |3| yields the measurement equation for the calibration setup: 



where we measure the distance R between the center of the transmitting antenna and the center of 
the lowest dipole of the Small Black Spider LPDA. 

4.4 Calibration Measurement Results 

In Fig. the absolute value of the amplified VEL of the Small Black Spider is displayed for the 
zenith direction d = 0° as a function of frequency. For frequencies lower than 30 MHz and higher 
than 80 MHz the reception is strongly suppressed due to the presence of the filter elements in the 
amplifier (cf. Sec. [OP . Within the bandwidth the VEL decreases with increasing frequency. This is 
the expected behavior for antennas which feature a gain that is constant as a function of frequency, 




(4.5) 




(4.6) 




Figure 10: The VEL of the Small Black Spider LPDA as a function of zenith angle for three different 
frequencies in measurement and simulation. Data points are given for all zenith directions that we 
were able to access with the setup shown in Fig. || . The uncertainty of the measurement is indicated 
by the bars analog to Fig. ^. 

such as LPDAs (cf. Eq. |A.8| ). The additional variations of the VEL within the frequency band occur 
due to the interplay of the LPDA's dipole elements which resonate at different frequencies. We 
observe that the simulation reproduces the bandwidth and the overall size of the measured VEL. 

For a set of three frequencies, the dependence of the VEL on the zenith angle is shown in Fig. 
[TD| . For the low frequencies, the antenna is most sensitive to zenith angles around 45°. At higher 
frequencies, a side-lobe pattern evolves with up to two lobes at the highest frequencies. 

The primary cause for the side lobes is the constructive and destructive interference of the direct 
wave and the wave reflected from the ground at the position of the antenna with its lowest dipole at 
a height of 3 m. Note that a conclusion on the reception of transient signals can only be drawn if the 
wide band combination of these patterns including their respective phasing is regarded, as we will 
do in Sec. P. 

With respect to the shape of the side lobe pattern we note a good agreement between measure- 
ment and simulation. For the combination of all zenith directions and all frequencies within the 
bandwidth we observe an agreement of the simulated and measured VEL of better than ±20%. 

The phasing of % a $ reveals the group delay induced by the Small Black Spider to the transmit- 
ted signal. The group delay X is given by: 

T(ffl) = -^arg(H^) . (4.7) 

In Fig. [TT] the group delay for the zenith direction is displayed. Within the measurement bandwidth 
from 30 to 80 MHz the group delay decreases by ~ 50 ns, where measurement and simulation agree 
on the functional dependence. 

As will be discussed in detail in Sec. |5] a non-constant group delay induces the dispersion of 
the observed transient signal and thus reduces its peak amplitude. The group delay displayed in 
Fig. [TT] results from a combination of the delay introduced by the logarithmic periodic structure of 
the Small Black Spider and the delay introduced by the filter elements of the amplifier. Although 
the amplifier has been designed to suppress the signal reception outside the bandwidth, it causes a 
non constant group delay also within the measurement bandwidth especially at lower frequencies. 
The most recent version of the AERA readout electronics is able to compensate digitally for the 
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Figure 1 1 : The group delay including the LNA for the zenith direction 6 = 0° as a function of 
frequency for the Small Black Spider LPDA in simulation and measurement. The uncertainty of the 
group delay is < 1.3 ns and results from variations observed in multiple measurements. 




Figure 12: The group delay of the Small Black Spider LPDA as a function of zenith angle for three 
different frequencies in the measurements and the simulations. 



dispersion induced by non constant group delays of the readout chain which we will report in a later 
article. 

In Fig. |1^ the group delay for the Small Black Spider is shown as a function of zenith angle 
for three different frequencies. We also find that the group delay exhibits a side-lobe pattern that is 
similar to the pattern observed in the case of the absolute values of the VEL. The shape of the group 
delay pattern is similar in the measurements and the simulations. The absolute values differ by up to 
30% depending on the considered frequency range, as is also visible in Fig. |TT]. 

Correction for a group delay that changes with incoming direction is important as single 
transient signal fronts will be observed at different detector stations. 

In the present calibration campaign we have measured the e$ -component of the VEL. The full 
VEL will be obtained in later calibration measurements. 

From the comparison of the measured and the simulated antenna characteristics we conclude that 
overall the simulations give a realistic description of the Small Black Spider antenna characteristics. 
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Figure 13: The simulated vector effective length as a function of frequency for the zenith direction 
for the three tested antennas. Note the logarithmic scale on the y-axis. 



5. Comparison of Transient Antenna Responses 

In this section we use simulations to compare the three antennas presented in Sec. with respect to 
their response to transient signals. 

5.1 Simulation of the Vector Effective Length 

To investigate the transient response characteristics of the antennas we simulated the respective 
VELs with the NEC-2 program. In the case of the Small Black Spider, the simulated antenna 
model is identical to the one presented in Sec. |]. Simulation models for the Butterfly and the Salla 
antennas have been created following their corresponding structure specifications. The Butterfly 
antenna is explicitly designed to be used at a height of 1 .5 m above ground, which we have set in 
the simulation. The other two antennas are simulated with the lowest sensitive element at 3 m, the 
current installation height of the AERA antennas. 

For the simulated calibration measurements in Sec. |] we have explicitly introduced a transmit- 
ting antenna as the source of the electric field. For the comparison between the different antennas we 
access the VEL from direct simulations of the receiving antenna without an additional transmitter. 
The simulated antenna characteristics presented in this section thus avoid possible deficiencies of 
the transmitter simulation and would correspond to calibration measurements performed at large 
distances. 

We use the antenna simulation program NEC-2 to access the VEL directly by exciting the tested 
antenna with a plane wave. Here the simulation computes the resulting currents in the receiving 
antenna which enables calculations of the response voltage as described in Sec. For each 
incoming direction, the two components of the VEL are accessed as the ratio of the response voltage 
and electric field amplitude by choosing the polarization of the plane wave accordingly. 

Alternatively, the tested antenna can be used as transmitter to access the VEL. With the NEC-2 
program it is possible to directly compute the far field generated by an antenna at a certain distance 
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Figure 14: The simulated group delay induced by the antennas for signals coming from the zenith 
direction 6 = 0°. A deviation from a non-constant group delay induces dispersion of the measured 
transient signal. 



R. By analogy to Eq. |3] the vectorial electric field emitted by an antenna is related to the VEL as: 

1 



£(R) = -iZ Q —l' Q He 



-imR/c 



(5.1) 



where Z( } is the current used to operate the antenna. The reciprocity theorem states the equivalence 
of antenna characteristics obtained from receiving and transmitting measurements. Solving Eq. |5TT] 
for the VEL yields the desired sensitivity of the antenna to incoming signals. Here it should be noted 
that the electric field £' given by the NEC-2 simulation is normalized to a unit distance of r = 1 m: 



£' = 8{R)Re iaR l c 



(5.2) 



We find that both options yield equivalent results for the VEL %. Using Eqs. 3.18 and 3.20| we 
finally produce % a which takes the impact of the amplifier and intermediate transmission lines into 
account. 



5.2 Characteristics of the Ultra- Wideband Vector Effective Length 

In Fig. [D] the absolute values of the VELs and in Fig. [B| the corresponding group delays are 
displayed as a function of frequency for the three tested antennas. In the case of the Small Black 
Spider, the depicted characteristics correspond to the displays in Figs. |9| and [TT] with slight changes 
due to the adapted simulation setup. 

The absolute values of the VEL are rather different for the three antennas. The amplifier of the 
Butterfly antenna exhibits the strongest gain resulting in large values for the VEL. The amplifiers of 
the Salla and the Small Black Spider feature similar amplifications. However, the VEL of the Salla 
is reduced as the antenna is loaded with an ohmic resistor. Note that the average levels of the VELs 
do not necessarily reflect the signal to noise ratio obtained in measurements. 
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Regarding the dependence of the VEL on frequency, both the Butterfly and the Small Black 
Spider focus on the performance at lower frequencies. The two antennas have been optimized with 
respect to the frequency content of the radio pulse which is predicted to be governed by longer 
wavelengths especially when observing air showers at large distances from the shower axis [43,44]. 

In addition, the Small Black Spider acts as a sharp bandpass to the AERA band. Also, the 
Butterfly slightly attenuates signals below 30 MHz to avoid short band transmitters. 

The drop-off in sensitivity visible for the Butterfly antenna at the highest frequencies results 
from the specific observation direction of 8 = 0° chosen for the display. Here, the constructive 
interference of the direct and the wave reflected from the ground is diminished at the height of the 
Butterfly of 1.5 m above ground for frequencies above 70 MHz. 

The Salla antenna has been designed for highest sensitivity at 70 MHz by approaching a 
situation of conjugate matching between the impedance of the antenna structure and the input 
impedance of the LNA (cf. Appendix 0). As will be discussed in Sec. the intensity of the radio 
background noise decreases with frequency. In combination with an increasing antenna sensitivity, 
the Salla is intended to measure noise spectra that are flat as a function of frequency. 

In Fig. ||5 the group delays of the antennas are pictured. In comparison to the Small Black 
Spider, the Butterfly antenna and the Salla feature an almost constant group delay within the 
considered bandwidth from 30 to 80 MHz. 

5.3 Transient Antenna Characteristics 

Up to now we have discussed the Fourier transform of the VEL H a . Antenna characteristics 
inspected at discrete frequencies correspond to the properties of the antenna that will be observed 
under the reception of mono-frequency signals which implies an infinite signal duration. The 
transient response is encoded in the development of the antenna characteristics as a function of 
frequency. 

To evaluate the distortion introduced to transient signals by variations of the VEL within the 
bandwidth we apply the inverse Fourier transform and inspect the VEL in the time domain: 

H a {t) = T-\U a {(0)) [mHz] . (5.3) 

Examination of the VEL in the time domain combines the full frequency range and takes into 
account the respective phase relations. Since the bandwidth is limited by the antenna characteristics 
the treatment of the VEL in the time domain results in a useful estimator. 

The result of the transformation in Eq. |53] is displayed in Fig. |T3| restricted to a common 
bandwidth of 30 to 80 MHz. In the time domain the antenna characteristics can be interpreted in 
terms of a transient wave form. The response to an incoming signal is calculated as the convolution 
of the respective displayed VEL and the incoming wave (cf. Sec. ^TTJ). 

For comparison, a bandwidth limited Dirac pulse is also shown. The shape of the Dirac pulse 
corresponds to a flat transfer function without dispersion. Hence, the Dirac pulse represents the 
VEL of an antenna that fully reproduces the shape of incoming signals within the bandwidth. 

The Salla and the Butterfly antenna have an almost unbiased response to incoming signals. 
The relatively large average group delay of the Small Black Spider, already indicated in Fig. [14], is 
visible in Fig. [T5] as a time delay of H a ^(t). In comparison to the other antennas the Small Black 
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Figure 15: The VEL in the time domain for the zenith direction 0=0°. The VELs have been 
limited to a common bandwidth of 30 to 80 MHz. The functions are renormalized to their respective 
maximum peak value which would be realized without dispersion in the antenna. The bandwidth 
limited Dirac pulse corresponds to an ideal antenna which introduces a bandwidth limitation to 
incoming signals only. 



Spider introduces the largest signal distortion. The pulse is broadened which corresponds to the 
dispersion induced by the variations of the group delay within the bandwidth. 

As normalization to the characteristics displayed in Fig. |T3| we choose the maximum absolute 
value of the respective H a $ (?) that would be realized if the VEL was not distorted by variations 
of its group delay. With respect to these maximum values we find that the Salla and the Butterfly 
antenna keep > 95% of their peak amplitudes whereas the maximum VEL of the Small Black Spider 
is reduced to ~ 70 % of the undispersed waveform. 

In the calibration measurements with the Small Black Spider we have seen that the directional 
properties of an antenna depend on the frequency under consideration (cf. Fig. |T0|). Extending 
the idea in Ref. [45], we investigate the maximum peak of the absolute value of the VEL in the 
time domain as a function of incoming direction. The result is displayed in Fig. [T^. The plot is 
subdivided for the two components of the VEL H a g and H a A. In the case of the Small Black Spider, 
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Figure 16: Peak directional diagram: The development of the maximum absolute value of the VEL 
in the time domain (cf. Fig. [T3J) is displayed as a function of zenith angle. The two components of 
the VEL are treated separately. The zenith angle dependence of the eg -component of the VEL is 
shown in the left part of the diagram. In the right part the -component is depicted. Two azimuthal 
directions have been chosen to maximize the readings of the respective VEL component (cf. Fig. 

The readings are normalized to the maximum peak value obtained in the zenith angle range. 
For the zenith direction 0=0° both components represent the same antenna characteristic which 
results in a connecting condition for the left and the right part of the diagram. 

it follows that the function H a $ (?peak> = 270°) in the right part of the diagram is a wideband 
representation of the directional characteristics displayed in Figs. [R] and [T2[ 

For the Small Black Spider antenna we find that the side lobe structure is less distinct in the 
case of the transient characteristic when compared to the side lobes at single frequencies (Fig. ITDI ). 
The peak response of the Salla antenna is reduced for zenith angles towards 6 ~ 50° when the 
e$ -component is considered. In this polarization direction the Butterfly antenna features the least 
complex coverage of the zenith angle range up to 70°. The suppression of zenith angles towards the 
horizon depends on the heights of the antenna above the ground. 

In the discussion of Fig. |5] it was emphasized that dipole-like antennas become insensitive when 
the incoming direction of the signal aligns with the dipole axis of the antenna. The development of 
the eg -component of the peak response is dominated by this geometric effect in the case of the Small 
Black Spider and the Butterfly antenna, as illustrated in Fig. [T5|. Here, the circular construction of 
the Salla antenna leads to a rather constant peak response up to high zenith angles. This results in an 
enhanced reception of vertically polarized signals from directions close to the horizon. 

From the simulation of the antennas we conclude that the Salla and the Butterfly antenna 
are best suited for the detection of transient signals as they introduce the fewest distortions to the 
recorded waveform. 
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Figure 17: Continuous spectra recorded simultaneously with the three tested antenna at the test 
bench of Nancay Radio Observatory over a 4 days period. An offset value is chosen for each 
antenna to allow for a common scale of the power axis. From top to bottom: Butterfly P ffset = -50.1 
dBm/MHz, Small Black Spider P offset = -62.9 dBm/MHz, Salla P offset = -67.5 dBm/MHz. See text 
for details. 

6. Comparison of the Reception of the Galactic Noise Intensity 

The response voltage to an incoming transient signal needs to be discriminated from a continuous 
noise floor. The sensitivity of an antenna can therefore be expressed in terms of signal to noise 
ratio as obtained from measurements. We have investigated the contribution of the signal to this 
ratio in the previous sections. To characterize the noise that is recorded with the antennas we have 
performed dedicated test measurements which are presented in this section. 

Beyond human-made noise the dominant source of continuous radio signal in the bandwidth 
from 30 to 80 MHz is the galactic radio background. This diffuse galactic background emission is 
mainly caused by charged particle gyrating in the magnetic fields of our galaxy [46]. 

The spectral irradiance / v (v) of galactic background depends on the frequency and the observed 
direction in the sky. The frequency dependence can be described by a power law: 

/ v (v)«v~P , [7]=Wm- 2 Hz- 1 . (6.1) 

The spectral index j8 ranges from 2.4 to 2.9 and depends on the frequency range and on the 
observation direction [47]. 

The dominant contributions to the galactic radio background come from the galactic plane. The 
tested antennas simultaneously observe a half-sphere of the sky. Consequently, variations of the 
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recorded noise level will occur due to the rise and fall of the galactic plane in the field of view of 
the antennas. The galactic noise background represents an irreducible background contribution to 
the observation of air shower signals. However, the variations of the galactic radio signal can be 
considered as a measure of the antenna sensitivity. These measurements were performed at the 
Nancay Radio Observatory. 

6.1 Observation of the Galactic Radio Background 

The antenna test bench at the site of the Nancay Radio Observatory is dedicated to the comparison 
and improvement of antennas used in radio astronomy. Owing to its radio quietness, the Nancay site 
features similar observing conditions as the site of the Pierre Auger Observatory. 

The three antennas under consideration (cf. Sec. 0) were installed in the test bench for 
simultaneous observation of the radio background intensity with their read out dipole axes pointing 
in east-west direction. Besides different coaxial cables, all antennas were read out with the same 
chain of analog electronics. Frequency spectra from 22 to 82 MHz were recorded from the different 
antennas in consecutive sweeps with the same spectrum analyzer. We measure one spectrum every 
~ 4 seconds per antenna. Within the measured bandwidth the continuous noise floor of the readout 
chain is more than 10 dB lower than the average signal power delivered by the Salla antenna, which 
provides the lowest signal power of the three antennas. 

In Fig. |l^the dynamic spectra derived from 10 5 spectrum sweeps for each antenna are displayed. 
The readings of the spectrum analyzer are normalized to the resolution bandwidth of the device to 
yield the recorded power spectral density. 

In the dynamic spectra, FM band transmitters appear as horizontal lines at the highest frequen- 
cies. At the lowest frequencies, short wave transmitters are visible with varying amplitude due to the 
changing conditions of the ionosphere. In the relevant frequency range from 30 to 80 MHz only a 
few transmitters are present, some of them only appearing for short periods of time. Vertical glitches 
are indicators for a saturation of the antenna over short periods of time. 

For all three tested antennas, a variation over a wide frequency range with a day-like periodicity 
is visible. This variation results from the changing position of the galactic disc in the field of view 
of the antennas and follows sidereal time. The sidereal day is ~ 4 minutes shorter than the Julian 
day. We have verified the observation of the sidereal period with repeated measurements in the test 
bench several months apart to rule out man-made interference [19]. 

We observe that the variation of the visible galactic radio sky is mapped differently by the 
antennas onto the recorded power spectral density. Especially in the case of the Salla, the amplitude 
of the variation appears to be less distinct than for the other antennas. In Sec. we will use 
the amplitudes to determine the fractions of noise that are added by the antennas internally to the 
recorded power spectral densities. 

For the site of the Pierre Auger Observatory, long term background observations are described 
inRef. [48]. 

6.2 Simulation of Galactic Radio Background Reception 

For proper judgment of the measured variations in Fig. [17], the directional sensitivities of the 
antennas need to be taken into account. We calculate a prediction of the progression of the power 



-26- 



J2000 




60 120 180 240 300 360 
Right ascension a [deg] 



Figure 18: Map of galactic noise intensity generated with LFmap at 55 MHz. Temperatures have 
been translated to intensities following the Rayleigh- Jeans law. The colored data show noise 
intensities in the field of view of the Nancay Radio Observatory at a specific time. The horizon 
is displayed as a dashed line, the star symbol denotes the local zenith direction. The shade over 
the colored data indicates the relative antenna sensitivity of the Small Black Spider LPDA at the 
corresponding frequency oriented in east-west direction at Nancay. The measured side lobe in Fig. 
[T0| (middle) is here pointing in the direction of the galactic center (Sgr A*). Please refer to a colored 
version of this plot. 



received due to the galactic radio emission using the LFmap program [49]. LFmap combines radio 
background maps measured at different wavelength and allows for an interpolation to arbitrary 
frequencies using appropriate power law indices (cf. Eq. |BTl"D . We obtain the background brightness 
B as maps of brightness temperature Tg in equatorial coordinates. The brightness is the spectral 
irradiance per unit solid angle and is given by the Rayleigh- Jeans law: 

2k 

B(v,a,S) = -^v 2 T B [S] =Wm- 2 sr- 1 Hz- 1 , (6.2) 
c z 

where k is the Boltzmann constant. In comparison to alternative descriptions of the galactic radio 
background given by Cane [50] and the global sky model GSM [51] we find an agreement of the 
three models at a level of ~ 1 dB by considering the spectral irradiance integrated over the full sky. 

An exemplary map of the galactic radio background generated at 55 MHz is displayed in Fig. 
[H| Besides the brightness, the dashed curve indicates the field of view that contributes to the 
recorded noise power at the location of the Nancay Radio Observatory at the given time. The relative 
directionality of the Small Black Spider antenna for the corresponding frequency is represented by 
the gray shade. 

The noise that is recorded in the measurement is the convolution of the currently visible radio 
background and the projection of the antenna characteristics onto the sky. The received power 
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Figure 19: Galactic Variation at the Nancay Radio Observatory at 55 MHz depicted for the Butterfly 
(top), Small Black Spider (middle) and Salla (bottom). The symbols display the measured data, the 
solid lines show the corresponding simulations. 

spectral density is given by: 

P v (v) = ^[B{v,a,8)A e (v,e,(j))da , (6.3) 

where A e is the effective aperture of the antenna including mismatch effects and losses (cf. Eq. |A.2| ). 
The factor 1/2 is introduced explicitly in Eq. to take the polarization mismatch between the 
antenna and an unpolarized radio brightness into account [52]. Using the result from Appendix |F], 
the recorded power spectral density can also be expressed using the VEL: 

Mv) = ~jB(v,a 1 S)\'H a (v,d,4))\ 2 dQ. . (6.4) 
a 

We have compared the measured and simulated average power spectral density as function of 
frequency based on Eq. using the same simulated VELs as for the discussion of the transient 
response in Sec. |5[ We find that the simulations reproduce the shape of the spectrum. However, 
they underestimate the noise intensity measured in the setup at the Nancay Radio Observatory by 
~ 4 dB for all three antennas. For the Small Black Spider we have cross checked the results of the 
simulation with spectrum observations at the site of AERA. Here, we do not observe differences of 
the absolute scale of measurements and simulations. In the following section we focus on relative 
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Figure 20: The variation of the galactic noise level obtained from simulations for the three antennas. 
The Butterfly is predicted to observe the largest difference in received power between galactic 
maximum and minimum. The variation of the Small Black Spider depends most strongly on the 
frequency. The galactic noise measured with the Salla is predicted to be least dependent on the time 
of observation. 



variations of the background noise as a function of sidereal time. Hence, the absolute values of the 
simulated noise are irrelevant for our further discussion. 

6.3 Comparison of Radio Background Variation 

In Fig. the variation of the power spectral density as a function of the sidereal time is shown for 
the measurements and the simulations. 

The measured variation corresponds to a cut in the dynamic spectra of Fig. [FT] at 55 MHz where 
the data have been condensed in the given time binning. For comparison, the measured and the 
simulated curves are shifted to a common mean value. 

The simulation overestimates the sidereal variation of the noise level without any additional 
noise sources but the galactic noise. The measured variation of the Butterfly antenna almost follows 
this ideal curve. If we assume that the simulated variation yields a valid prediction of the noise 
progression, we conclude that additional external noise sources are negligible for our measurement. 
At 55 MHz the Butterfly realizes a variation of ~ 2.5 dB between galactic maximum and minimum 
where the simulation predicts a maximum variation of ~ 3.1 dB. For the other two antennas the 
variation is smaller. 

To evaluate the measured variations with respect to noise that is added in the antenna internally 
we have to take into account that the predicted maximal variation is different for each antenna. 
Displayed in Fig. ^ is the maximum variation for the full frequency range. 

The projection of the directional properties of the Butterfly antenna yields a larger maximal 
variation than in the case of the other two antennas. The deviations throughout the bandwidth in the 
case of the Small Black Spider indicate the lobe structure of the antenna evolving with frequency. 
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Figure 21: The measured galactic variation for the three antennas as a function of frequency here 
indicated with the bar symbols. The continuous lines give the expected variation from simulations 
as in Fig. ^0] but with a noise power added that is constant in frequency. The percentages denote 
the strength of the added noise floor relative to the maximum signal expected from the galactic 
background. 



The displayed maximum variation V m ax is the ratio of maximal to minimal received galactic 
noise: 

W(v) = j-+ ■ (6.5) 

For simplicity we assume that the antenna is adding internally an additional noise floor Pj nt that is 
constant with frequency. The realized variation is then: 

Vr(v) = W(v)+P mt 



We are interested in the fraction / of internal relative to the galactic noise. In Eq. q]6J we count the 
fraction / relative to the maximum noise that is predicted throughout the whole bandwidth: 

v i \ _ Anax ( v) + / ' -Pmax ( Vmax) 

I{ '"Uvl+Z'Wv.-.) ' 

In Fig. ^T] the measured variation is displayed as a function of frequency for the three antennas. The 
simulated variations from Fig. ^3 have been adjusted with noise fractions / to match the realized 
variations. 

The observed variation of the Salla antenna ranges below 1 dB. Therefore, even at the galactic 
maximum, more than half of the output power corresponds to internal antenna noise. 

The Small Black Spider almost realizes the full variation at the lower frequencies indicating a 
fraction of internal noise of ~ 25%. However, at the highest frequencies the variation of the galactic 
background is less distinct. Here, the antenna's sensitivity is not sufficient to raise the galactic signal 
strongly above the internal noise floor. 
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The Butterfly antenna provides a power spectral density which is strongly dominated by 
the galactic noise background over the full frequency range. This constitutes the best observing 
conditions to air shower signals in terms of continuous background noise. 

7. Summary and Conclusions 

With AERA — the Auger Engineering Radio Array — the Pierre Auger Collaboration addresses 
both technological and scientific questions of the radio detection of ultra-high energy cosmic rays. 

AERA will check the feasibility of the radio-detection technique on large instrumented sites 
and is thus intended to serve as a blueprint for the next generation of ground-based cosmic ray 
detectors at very high energies. 

The installation of the first stage of AERA (21 stations) is the startup for the construction of 
a 20 km 2 radio detector consisting of 160 autonomous self-triggered detector stations sensitive to 
frequencies from 30 to 80 MHz. The first stage has been taking data since September 2010. 

Along with the construction of the first stage of AERA, we have performed dedicated studies 
of the antenna acting as sensor to the radio emission. The transient nature of the air shower 
signal requires a detailed description of the antenna sensor, aiming for a calibrated measurement 
of the incident signal. We identify the vector effective length as suitable quantity to perform the 
calculation of the antenna response to transient signals including multiple reflections, interference 
and polarization effects. 

Having identified the relevant antenna characteristics we reported the calibration of a logarithmic- 
periodic dipole antenna used at AERA. For the first time, the zenith angle dependency of the sensi- 
tivity of a full scale radio detector antenna was accessed including the vectorial phase information. 
This allowed us to calculate the vector effective length of the antenna on an absolute scale. 

For the next setup stage of AERA we have evaluated three candidate antennas each pursuing 
a different strategy for an optimal reception of cosmic ray signals. A logarithmic-periodic dipole 
antenna, a loaded dipole loop antenna and an active bowtie antenna were considered. 

On the basis of the vector effective length we studied the transient response characteristics of 
the antennas in the time domain. Our analysis reveals a considerable reduction of the peak amplitude 
response due to dispersion in the case of the logarithmic-periodic dipole antenna. 

At the Nancay Radio Observatory we have performed simultaneous measurements of the 
galactic radio background with the three antennas. Observing the variation of the measured noise 
power as a function of sidereal time we have tested the antennas for an optimal continuous noise 
performance. 

We find that the active bowtie antenna is superior in both the noise performance and the transient 
response. For the next setup stage of AERA we have chosen the bowtie as sensor to the coherent 
radio emission from cosmic ray induced air showers. 
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A. Vector Effective Length and Gain 

According to Eq. the VEL allows antenna calculations to be performed on the basis of complex 
amplitudes. The VEL is related to more commonly used quantities to describe antennas such as the 
antenna directivity. Directivity is related to the induced signal power. In this section we summarize 
the interrelation between the power- and the amplitude-based calculations. 

The intensity of a wave averaged over time in a single polarization direction St with k = 6 , <j> at 
a specific frequency is given by: 

ft-i^e (a.i, 

2 Z 

with Zo ~ 120 % Q. the vacuum impedance. The power that is available to the readout (load, L) of 
the antenna is then accessible via the maximum effective aperture A em ,k of the antenna: 



= Aemi 



S k ■ (A.2) 
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Figure 22: Representation of the vector field of the VEL depicted in Fig. g as the antenna gain. For 
each incoming direction of a signal the gain is given as the distance between the surface and the 
center of the coordinate system on a logarithmic scale. Left: The components of % in ee-direction 
yield the partial gain Gg which is sometimes called vertical gain. The dashed blue curve presents a 
cut in the gain sphere which is referred to as E-field plane gain. Right: The components of % in 
e<Sp -direction yield the partial gain G$ also referred to as horizontal gain. Here, the blue dashed curve 
is called H-field plane gain. 

In antenna theory [38] the effective aperture is related to the directivity D of the antenna. The 
directivity for a specific polarization is called partial directivity: 

A 2 

A Qm ,k = ~rD k . (A.3) 
471 

Hence the power that is delivered by the antenna structure is: 

A 2 I £k I ^ 

/, «=8i D ^ • (A " 4) 
Using the partial response voltages from Sec. [O] we introduce the VEL into the power calculation: 

. (A,, 

87T Zo 

The directivity describes the relative antenna properties to different incoming directions of the signal 
only. It does not take into account losses inside the antenna structure nor possible mismatches 
between the antenna and the readout system. The impact of such mismatches is discussed in Sec. 

The power given in Eq. |A.5| will only be achieved under the condition of conjugate matching. 
Conjugate matching is discussed in detail in antenna literature as it describes an optimal condition 
for the reception of signals and simplifies the relevant equations [38]. In this case the power available 
at the readout is related to the open circuit voltage V c that is induced by the signal over the footpoint 



impedance Za of the antenna: 



V -' 2 



rL ' k= SRe(z7) ■ (A " 6) 
The standards in Ref. [53] define the VEL to map the electric field to the open circuit voltage. Hence 
we identify V c,/t with the partial voltages Ve and mentioned in Sec. |3~T| . The combination of Eq. 
|A.5| and Eq. |A.6| then yields the relationship between VEL and the directivity: 

d ^r^T^ ■ (A7) 

When accessing the VEL in simulations and measurements it is usual to include losses inside 
the antenna structure e.g. due to ohmic resistance. In terms of power, these losses are accounted for 
by multiplying the directivity with a dimensionless efficiency factor e. If losses inside the antenna 
structure are included in the VEL, Eq. |A.7| relates to the antenna gain G: 

G <= eD '=«lb£ |Wi|2 ■ <A - 8) 

The treatment of losses inside the antenna structure is important as the loaded loop antenna Salla (cf. 
Sec. 0) explicitly introduces an ohmic resistor into its structure to give a specific shape to its gain. 

It should be noted that here the gain is a real number for each of the two polarization directions. 
Hence the gain is not sufficient to describe subtle effects such as the polarization dependence of the 
recorded signal as discussed in Sec. or dispersion effects within the antenna structure which 
invoke a complex phase. To save the concept of a power based calculation, a phase can be introduced 
in addition to the gain and polarization mismatch factors can be added to the calculation. 

In Fig. ^ the partial gains corresponding to the T~Lq and the Ti^ component of the VEL at a 
single frequency are displayed. Consequently, the two plots present a subset of the information 
given in Fig. For a specific direction (0,0) the magnitude of the gain is displayed as the distance 
from the center of the coordinate system. The gain spheres given in Fig. and subsequent cuts 
through the patterns are more common representations of antenna characteristics. They are helpful 
when discussing the directional properties of antennas. 



B. Reconstruction of the Electric Field in Dual Polarized Measurements 

In current radio detection experiments, there are at least two antennas installed at each observing 
position. In this section we will discuss a major benefit of such dual measurements. That is the 
reconstruction the 3 -dimensional electric field vector of the incoming signal. 

In the case of the antenna models discussed in Sec. 0, two rotated antennas are assembled in 
the same hardware structure (cf. Figs. 0, [3] and 0). Typically, one antenna is rotated by 90° with 
respect to the other. In principle, such a setup allows for two independent measurements of the same 
electric field 6 with two different VELs Hi and H2. The corresponding vector field of the combined 
measurement is depicted in Fig. ^3|. It results in two recorded signals that relate as follows: 

Vi(o) =Hi (<M,0) •£(<») , 

= fti, e (a>,e,0)£fl(fi>) +n^{(o,6,^)£ ((o) (B.l) 
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Figure 23: The VEL in the case of a dual polarized antenna setup. Due to two independent 
measurements of the same incident electric field, the full vectorial signal information can be 
reconstructed when the incoming direction of the signal is known. 



and 

V 2 (o) = H 2 (<o,e, <!>)■£ {co) , 

= U 2 ,e{(D,d^)£ e {(D) +n 2 .^co,d^)£^(o) . (B.2) 

The system of Eqs. |B.1| - |B.2| can be solved for the two components of the electric field vector. To do 
so, the antenna characteristics need to be evaluated at the incoming direction (0, 0) of the signal. 

Here it should be noted that the wave front emitted by an air shower deviates from a planar 
wave front [54]. Hence, the incoming direction of an air shower signal at an antenna sensor is 
expected to slightly deviate from the direction of the air shower axis. In Ref. [55] we find that this 
introduces only minor effects in the reconstruction. Changes in the antenna sensitivity are typically 
small on the scale of angular variations of the incoming direction, which range at a few degree. 

Using, for instance, the direction of the air shower axis the system of Eqs. |B.1| - |B.2| can be 
solved yielding the full electric field vector: 

6{ ' n l6 (<o)n 2 ^(w)-n^((o)n 2 , e ((o) K ' } 

V 2 (co)-n 2 , e {co)£e(co) 

£ * {(o) = ?M5) • (B - 4) 

An implementation of this reconstruction scheme is implemented in the 'Offline' software framework 
of the Pierre Auger Observatory [56]. 

At this point the response of the dual polarized measurement setup can be expressed in terms of 
a matrix calculus when the two measured voltage functions in Eq. |B.1| are interpreted in terms of a 
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response vector: 

V(co) = M .£(co) . (B.5) 

For single frequencies the matrix containing the components of the VEL is referred to as Jones 
antenna matrix or voltage beam matrix in radio polarimetry [31]. 



C. Multiple Reflections in Intermediate Transmission Lines 

A circuit diagram for a setup with intermediate transmission line is displayed in Fig. [7] (right). Here 
the open circuit voltage induces an initial forward traveling voltage V + into the transmission line of 
impedance Z t \. Using the result from Eq. [3.17| this voltage is: 

V + = -4^Voc • (C.l) 

Looking at Fig. we derive a model for the desired voltage amplitude over the readout impedance 
Vl- The complex voltage amplitude at the end of the transmission is the initial voltage V + multiplied 
with a propagation factor: 

V + (at/ / ) = V + -e- (! '? +a)/ ' = V + -e- 1 "' , (C.2) 

where y is the complex propagation constant and a the attenuation loss per unit electrical length /' 
of the transmission line. According to the Fresnel coefficients, the voltage over the load impedance 
follows as: 

V<u = V + -e-i"'(l+r L ) . (C.3) 



Here Yi is the voltage reflection coefficient: 



z L +z tl 



(C.4) 



Following Fig. 0, part of the amplitude is reflected back towards the antenna. There a further 
reflection 1^ occurs at the antenna impedance including the transformer operating with the inverted 
transformation ratio: 

Zi/r — Zfi 
Z A /r + Z tI 

If the scheme in Fig. ^ i s continued it yields the total voltage over the load impedance Vl as the 
sum over all amplitudes transferred from the antenna to the load at the end of the transmission line: 

DO 

V L = £ Vi, L (C.6) 

n=0 

= ...=v+(i+rj- ?37 ^ . (c.7) 

In Fig. ^5] the result for the response Vi(t) = J 7-1 (Vl(w)) to an exemplary open circuit voltage V c 
in the time domain is displayed. In this example a setup of Z& = Zl = 200 Q., Z t i = 50 £2 and an 
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Za/t 



v +e -^'(i + r L ) A 



V + e-^ 3l 'T L r A (l + T L ) 



V+ e-^ l 'T 2 L T 2 A (l + T L ) 



Figure 24: Model for multiple reflection between mismatched source and load impedance with 
intermediate transmission line. 



electrical line length of /' = 9 m is chosen. These parameters lead to a reappearance of the pulse at 
the load due to multiple reflections. As a cross check the result for the corresponding setup obtained 
with the circuit simulator program QUCS [57] is displayed showing excellent agreement. 

In summary, combining Eq. |3.15| and Eq. |C.7| we obtain the transfer function that will be 
realized in a measurement setup with the considered antennas by: 



VrZ tl 

Z A + rZ n -T A T L 



(C.8) 



D. Renormalization of S-Parameter S21 

When characterizing an amplifier with scattering parameters, the S-Parameter S21 yields the 
amplified voltage amplitude normalized to the incoming voltage amplitude V g delivered by the 
signal generator: 

S21 = -^ . (D.l) 

The realized VEL discussed in section ^3] relates the incoming electric field to the voltage over the 
input impedance of the amplifier Vl- From the S-parameter measurement the voltage at the amplifier 
input is calculated analog to Eq. |C.3| using the voltage reflection coefficient Fl- 

V L = V,(l+r t ) = V,(l+Sll) (D.2) 

where the reflection coefficient corresponds to the scattering parameter S 1 1. Hence, the amplification 
of the response of the antenna to an electric field calculated with the VEL is given as: 

V a V a S21 



V L V g (l + Sll) 1+S11 



S21' . (D.3) 
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Figure 25 : Analytic calculation of the voltage Vi over the load impedance as a response to an initial 
voltage pulse V oc - for the circuit diagram depicted in Fig. [7] (right). The calculation is performed 
with the ansatz presented in Fig. The result is compared to the response simulated with the 
circuit simulator QUCS [57]. 



E. Vector Effective Length and Realized Gain 

The active component of the power that is consumed by an antenna used as transmitter is given by: 

P = ¥lL.Re(Z* a ) . (E.l) 

Zq is the current within the antenna structure analog to the receiving case depicted in Fig. [7]. Z' a 
is the impedance of the transmitter. The power that is delivered through the coaxial cable to the 
transmitter is: 

r* = ^r ' (E2) 

where the voltage amplitude V g is the denominator of the S -parameter S21. Due to the impedance 
mismatch between the transmitter antenna and the coaxial cable, only a fraction of the power P g is 
available to be consumed by the antenna: 

\? Q f.Re(Z t a ) = ^(l-\r'\ 2 ) , (E.3) 

with r f the voltage reflection coefficient of the transmitting antenna to the coaxial cable. We resort 
Eq. |E.3| and multiply with the squared VEL of the transmitter: 

|T °' 2| ?4I 2 = D ,iw (l-m 2 )|74l 2 • (E.4) 



|V g | 21 0' Re(Z' a )Z tl 
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In Eq. |A.8| we have derived the relationship between VEL and gain. Introducing this interrelation on 
the right hand side of Eq. |E.4| yields: 



IT' I 2 
l x ol \qjt 

IV |2 1 



(1 



~t\2- 



G' 



(E.5) 



cal 



which is the result claimed in Eq. 4.5. The calibration of the transmitter antenna GLj is given by the 
manufacturer including the reflection at the input of the antenna. 



F. Effective Aperture and Vector Effective Length 

In Eq. |A.2| the power available to the load due to an incident wave is given in terms of the maximum 
effective aperture of the antenna. The maximum effective aperture is obtained under conjugate 
matching and without losses in the antenna. If losses and mismatching effects are included we refer 
to the effective aperture A s as: 

Pl=A & S . (F.l) 
The available power is also accessible via the realized vector effective length: 

D 1 |V L | 2 1 \u,x\ 2 



Z L 2 Z L 



(F.2) 



The vector product in Eq. is the projection of the electric field vector onto the vector of the VEL. 
To disentangle the vectorial from the power calculation a polarization factor p is introduced [27] : 



\n r £\ 2 

P \H r \ 2 \£\ 2 

such that: 



(E3) 



1 Z L 2 V 
Using Eq. |A. 1| we introduce the intensity of the wave which is also used in Eq. jFTTj : 

PL = P~\n r \ 2 S . (F.5) 

Zl 

Hence we derive the relation between realized VEL and the effective aperture: 

A e = p^\U r \ 2 . (F.6) 
Zl 

For an unpolarized signal the electric field vector and the vector of the effective length will align 
randomly in time. In time averaged measurements the expectation value of the polarization factor is 
hence: 

<p>=l/2 . (F.7) 
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G. List of Acronyms 



AERA Auger Engineering Radio Array 

AM Amplitude Modulation 

ASIC Application-Specific Integrated Circuit 

CODALEMA COsmic ray Detection Array with Logarithmic ElectroMagnetic Antennas 

FM Frequency Modulation 

FPGA Field-Programmable Gate Array 

GPS Global Positioning System 

LMC Large Magellanic Cloud 

LNA Low Noise Amplifier 

LOFAR LOw Frequency ARray 

LOPES LOFARPrototypE Station 

LOPES -STAR LOPES Self Triggered Array of Radio detectors 

LPDA Logarithmic-Periodic Dipole Antenna 

NEC Numerical Electromagnetics Code 

SALLA Short Aperiodic Loaded Loop Antenna 

SMC Small Magellanic Cloud 

VEL Vector Effective Length 

QUCS Quite Universal Circuit Simulator 
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